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ABSTRACT 

Aims. We seek to constrain the surface composition of the Trans-Neptunian Object (90482) Orcus and its small satellite Vanth, as 
well as their mass and density. 

Methods. We acquired near-infrared spectra (1.4-2.4 fim) of (90482) Orcus and its companion Vanth using the adaptive-optics-fed 
integral-field spectrograph SINFONI mounted on Yepun/UT4 at the European Southern Observatory Very Large Telescope. We took 
advantage of a very favorable appulse (separation of only 4") between Orcus and the UCAC2 29643541 star (m R = 11.6) to use the 
adaptive optics mode of SINFONI, allowing both components to be spatially resolved and Vanth colors to be extracted independently 
from Orcus. 

Results. The spectrum of Orcus we obtain has the highest signal-to-noise ratio to date, and we confirm the presence of H2O ice in 
crystalline form, together with the presence of an absorption band at 2.2 fim. We set an upper limit of about 2% for the presence of 
methane, and 5% for ethane. Because the methane alone cannot account for the 2.2 /jm band, the presence of ammonia is suggested 
to the level of a couple of percent. The colors of Vanth are found slightly redder than those of Orcus, but the large measurement 
uncertainties forbid us from drawing conclusions on the origin of the pair (capture or co-formation). Finally, we reset the orbital phase 
of Vanth around Orcus, and confirm the orbital parameters derived by Brown et al. (2010, AJ 139). 

Key words. Kuiper belt objects: individual: (90482) Orcus; Methods: observational; Techniques: high angular resolution; Techniques: 
imaging spectroscopy 

1. Introduction 2. Near-infrared spectroscopy of Orcus has revealed a 



surface rich in water ice in crystalline form (Fornasier et al. 
Binaries in the Solar System are of high importance because they 2004; de Ber g fad^j] Eool [Truiillo et al. 2005; Barucci et al. 
provide the most direct and precise way to derive the mass of I200M iDeMeo et al. 2010). Moreover, Truiillo_etalJ <|2005), 



minor planets (see | Hilton || 2002|) . Combined with volume esti- Barucci et al.l(l2008bl) . and Sil dZOlfJ ( detected a weak 



mates, their densities can be calculated, providing hints on their band 2 2 that could be associated with either met hane 

composition and interior (e.g., | Merline et alj | 2002t | BrittetalJ (CH4) Qr ammonia (N h 3 ). The long term stability of all ices 

|2002J). Subsequently, they constrain the characteristics of the are affected by high energy photon bombardment (causing pho . 

most pristine material of the solar system, and further our un- todissociation and sput tering), micrometeorite impacts, radioac- 

derstanding of planetary system formation and dynamical evolu- tiye decayi and sublimation . Both methane and ammonia ^ ex - 

tion. In this valuable context, the Trans-Neptunian Binary (TNB) p£Cted tQ be destroyed by so i ar irradiation on short timescales 

Orcus/Vanth system is of particular interest for the following rea- fetrazzulla & Palumbol fl998t ICooner et al.l l200l ICottin et al.1 

sons - e j-. , 12003b . Ammonia's presence, if confirmed, would thus require 

1. With an estimated albedo of -27% (|Lim et ajj |2010|), an actiye procesfj tQ resupply the surface with ammoniai such as 

Orcus is among the brightest known Trans-Neptunian Objects impact gardening or> more favorab i y , cr y OV olcanism (the ammo- 

(TNOs); and has a diameter of about 850 kilometers. nia lowers the melting tem perature of water ice and h ence favors 

such mechanism as highlighted bv lCook eta l. 2007). Bodies in 



Sendoffprint requests^ to: Benott Ca^: braoiLcai^Sesajnt ^ the QUter ^ system ^ haye methane on meir surface have 

retained their atmospheres which would also have in 
plications for its discovery on the surface of Orcus. 
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Based or. observations collected at the European Southern retained their atmospheres which would also have i mportant im . 
Observatory Very Large Telescope (programs ID: 284.C-5044 & 



384.C-0877I. 
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3. Recent radiometric measureme nts from ESA He rschel 
(Key Program "TNOs are Cool!", see iMuller et alJl2Q09h have 
refine d the size estimate of Orcus to 850 ±90 km (iLim et al.l 
l2010h . The diameter estimate will potentially be improved from 
the stellar occultations expected for upcoming ye ars. Thus the 
impro vement of Vanth's orbit (upon the solution bv lBrown et al.l 
2010) will help determine the bulk density of Orcus. 

We present here new spectro-imaging data obtained in 2010 
that provide constraints on the composition of Orcus and the or- 
bit of Vanth. We describe in Section [2] the observations, list in 
Section|3]the data reduction and spectral extraction steps, present 
in Section|4]the analysis of the colours and spectra of Orcus and 
Vanth, and detail in Section [5] the orbit computation and stellar 
occultations prediction. 

2. Observations 

The brightness contrast (Amy ~ 2.6) and small apparent angu- 
lar separation (~0.2") between Orcus and its satellite Vanth re- 
quire the use of a high angular-resolution camera/spectrograph 
to spatially resolve the system. This means observations have to 
be co nducted in the visi ble from the Hubble Space Telescope 
(e.g., iBrown et"al .1 120101) . or in the near-infrared with ground- 
based telescopes equipped with adaptive optics (AO) modules. 
The latter is of high interest for cold objects like TNOs because 
many ices display strong absorption bands in the near-infrared 
(see lBarucci et al.ll2008aL 

However, because adaptive-optics systems require a bright 
(my < 15) reference source (a Natural Guide Star: NGS) to 
correct the incident wavefront from the deformations induced 
by the atmospheric turbulence, study of TNOs from the ground 
with AO is generally limited to the brightest objects {e.g., Pluto 
or Haumea). The extension of such studies to fainter targets is 
possible thanks to two techniques: a laser beam can be pro- 
jected into the atmosphere to create an artificial star of mag- 
nitude Mr ~ 13.4, called a Laser Guide Star (LGS). However, 
because the laser beam is deflected on its way up by the atmo- 
spheric turbulence, the LGS position moves on the plane of the 
sky in a random pattern (corresponding to low orders of the tur- 
bulence, called tip-tilt). Hence, a natural close-by star must be 
monitored to correct the wavefront from the motion of the LGS. 
Because the requirement on these reference stars (called Tip-Tilt 
Star: TTS) are less strict (angular distance and brightness) than 
for NGS, several TNOs have already been observed this way 
(e.p jBrown eTaT1l2006t iDumas et al.ll20llh . 

Another solution consists of computing close encounters 
(separation smaller than about 30") on the plane of the sky be- 
tween the object of intere st and a star suitable as a NGS (e.g., 
Bert hier & Marchisll200ll) . These events are called appulses. On 
2010 February 23 UT, Orcus had a particularly favorable appulse 
with the star UCAC2 29643541 (m R = 1 1.6) at an angular sepa- 
ration of only 4". We thus observed it in Service Mode (program 
ID: I284.C-50441 at the European Southern Observatory (ESO) 
Very Large Telescope (V LT) with the near- infrared integral- field 
spectr ograph SINFONI (lEisenhauer et alJ l2003t iBonnet et all 
2004). Observations were realized simultaneously in the atmo- 
spheric H and K bands ( 1 .45-2.45 fim) using the H+K grating of 
SINFONI, providing a spectral resolving power/? of about 1500. 
We used a plate scale of 50x100 mas/pixel, associated with a 
3" x 3" field of view. We alternated observations of Orcus and 
nearby sky in a jitter pattern to allow optimal sky correction, be- 
ing cautious to avoid the NGS (4"). 

Unfortunately, the AO module of SINFONI has not been de- 
signed to offer differential tracking (i.e., NGS fixed on the plane 



of the sky, field of view following a target with non-sidereal mo- 
tion). We thus had to set the duration of integrations as a com- 
promise between the slew of Orcus on the detector plane and the 
count level reached on Vanth (my ~ 21.6). We used individual 
exposures of 150s to theoretically^ achieve an average signal- 
to-noise ratio of 1 on Vanth over H band. In return, during a 
single exposure, Orcus moved by -0. 109" in right ascension and 
0.039" in declination, distorting its apparent shape, thus elon- 
gated along the SE-NW direction as visible in Fig.Q] 

Atmospheric conditions at the time of the observations were 
very good, with an average seeing of 0.8" and a coherence time 
ranging from 7 to 20 ms. Orcus was close to zenith during the 
observations with an airmass ranging from 1.05 to 1.4. This al- 
lowed the AO system to provide optimal correction, resulting in 
a spatial resolution close to the diffraction limit of the telescope 
(Orcus FWHM was of 85 x 100 mas in K-band). 

We also report here on some test observations of Orcus 
performed on 2010 March 13 UT at the ESO VLT (prog. ID: 
384.C-0877) in the so-called "Seeing Enhancer" mode. This 
mode consists of closing the AO loop on a LGS, but without 
providing any TTS, Orcus itself being too faint (my ~ 19.7) 
to be used as TT S (as opposed to targets like Haumea, see 
Dumas et al. 2011, for instance). Hence, only the higher orders 
of the atmospheric turbulence are corrected (i.e., there is no tip- 
tilt correction). The advantage in that mode was to perform dif- 
ferential tracking and to take longer exposures (600s). The in- 
strument settings and observing strategy were otherwise similar 
to those for February observations. 

Atmospheric conditions were worse during March observa- 
tions, with an average seeing of 0.9", and coherence time of 
about 3 ms. However, the quality of the correction provided by 
the AO in that mode is intrinsically lower than in the case of 
the appulse: Orcus FWHM was of 0.38" (still representing an 
improvement by a factor of =2 with respect to seeing-limited 
observations). Despite the shape of Orcus being elongated by its 
apparent displacement in February, the quality of the data was 
superior (with a shorter exposure time) to the March data where 
the spread of its light is directly related to the lower AO cor- 
rection achieved. This highlights the advantage of searching for 
favorable appulses for faint moving targets to use bright NGS as 
reference for the adaptive optics correction. 

3. Data reduction and analysis 

We used the SINFONI pipeline (Mod idiani et al.ll2007l) version 
2.0.5 to perform the basic data reduction: bad pixel removal, flat 
fielding correction, subtraction of the sky background from the 
jittered observations, and wavelength calibratio n with Xenon- 
Argon-Krypton lamps (see Guilb ert et alJl2009l for a complete 
description of the procedure on other faint TNOs). Default pa- 
rameters were used except for in the "jitter" recipe where we 
set the parameters scales. sky to true and density to 3 for opti- 
mal sky-background correction. This provided us with 27 and 9 
individual cubes (two spatial plus one spectral dimensions) of 
Orcus/Vanth for February and March observations, respectively. 

We then computed the average centroid position of Orcus for 
each individual observation by stacking the cubes along wave- 
length. We used this information to shift and add all the cubes 
into a single one for each date, corresponding to equivalent expo- 
sure times of 4050s and 5400s. We then re-aligned all the wave- 
length slices of the cube because the centroid position of Orcus 
was not constant with wavelength but rather presented a slow 

1 computation made using ESO Exposure Time Calculator 
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2010 Feb 23 I 2010 Mar 13 




Fig. 1. Two images of the Orcus-Vanth system, obtained by sum- 
ming of all individual observations and stacking the resulting 
cube along wavelength. Left: 2010 February 23 UT observations 
using the appulse NGS AO correction, representing a total inte- 
gration time of 4050s. Both Orcus and Vanth are easily separa- 
ble in the image. Right: 2010 March 13 UT observations using 
the LGS AO correction, without TTS reference (see text), corre- 
sponding to a total integration time of 5400s. The angular resolu- 
tion provided in that mode forbids the detection of Vanth, which 
flux, spread over many pixels, is hidden within the background 
noise. 



dr ift due to th e differential atmospheric refraction as described 
in lCarrvetal1d2010h . 

We then extracted the respective spectra of Orcus and 
Vanth by adjust ing (using MPFIT least-square method from 
Markwardt 2009), for each wavelength, a model I composed by 
a linear background and two Moffat functions describing both 
components: 

I(x, y) = T (x, y) + T v (x, y) + ax + by + c ( 1 ) 

where Ti are the two Moffat functions, representing Orcus (T ) 
and Vanth (?>), defined by: 

1 + ((x-x^-iy-yiyM 2 
T i {x,y)=f i . *' °> I (2) 

+ J(x-j[J)^ + (y-yp^2) 

where x, y are the frame spatial dimensions, fi is the peak level of 
each Moffat function centered on the coordinates {x c j7 y?). cr x , cr y 
are the Half-Width at Half Max (HWHM) along two perpendic- 
ular directions, making an angle 8 with the detector x direction, 
and a is the power law index of the Moffat functions. The final 
spectrum was cleaned for bad points using a 3 cr median smooth- 
ing procedure. 

The advantage of this method is to provide the spectra of 
both components as well as their relative astrometry. We discuss 
both points in following sections. 

4. Spectral Analysis 

4.1. The surface composition of Orcus 

Figure [2] compares o ur new spectrum of Orcus to that of 
Baru cci et al.l d2008bl) . The overall spectral shape reveals the 
presence of water ice, dominated by the crystalline form 
as already addressed in previous work dFornasier et alj [ 2004 ; 
de Bergh et al I l2005b iTruiillo et al J l2005b IBarucci et al.ll2 008b: 
Guil bert et alJl2009t Delsanti et alJBO 1 0t DeMeo et al.ll20 1 Oh . 

The ratio of these two spectra, shown in the bottom part of 
Fig. [2] reveals their similarity, although we note a difference in 




1.50 1.75 2.00 2.25 

Wavelength (/urn) 



Fig. 2. Upper panel : Spectra of Orcus in th e H and K bands from 
this work (top) and lBarucci et al.l d2008bl) (middle). The bottom 
spectrum is an average of the two. The spectra are normalized to 
1.0 at 1.75/im and are shifted by +0.65, and -0.4. The spec- 
trum taken in March 2010 was significantly noisier than the two 
shown here so we do not plot i t nor do we us e it in our analysis. 
Lower panel : The ratio of the IBarucci et al.l data and data from 
this work which shows little difference between the spectra apart 
from a small flux difference shortward of 1.65 jjm. 



the overall flux level (~10%) shortward of ~ 1 .65 ^im. This differ- 
ence does not appear to be related to variation of H2O (amount 
or grain size) because it is present shortward of 1 .5 fim. Potential 
explanations include instrumental effects or differences of the 
standard stars. 

We confirm the presence of a feature near 2.2 /jm with a band 
center located at 2.209 + 0.002 yum and a band depth of 9 ± 2% 
that previous works have att ributed to CH4, NH3, or NH^. We 
combined the Baru cci et al.l spectrum with ours to slightly in- 
crease the ov erall signal-to-noise r ati o. We did not i nclud e the 
spectra from iDelsanti et akl d2010t) or iDeMeo etaH d2010h be- 
cause the quality of those data were significantly lower. This av- 
erage spectrum is used for all of the analysis reported in this 
section. 

Current volatile retentio n models (e.g., ISchaller & BrowrJ 
l2007t iLevi & Podolakll2009h predict that CH 4 is not stable on 
Orcus' surface over its lifetime, however, Orcus' intermediate 
size among TNOs place it closer to the retention boundary than 
most other objects and provides us with an opportunity to test 
these models and perhaps make constraints on the assumptions 
therein. An important step in understanding Orcus' surface com- 
position is thus a search for weak bands hidden in the spectrum 
near the detection limit. While many species could potentially 
exist in small quantities on the surface, the lack of multiple 
strong bands makes their identification difficult. Here we focus 
on searching for methane bands in the spectrum primarily be- 
cause of the abundance of strong bands in the appropriate wave- 
length regime, but also because of the important implications a 
detection would have on our understanding of the surface condi- 
tions of these bodies and the criteria for volatile retention. 
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Fig. 3. Plotted here is the average spectrum of Orcus from 
our work a nd Baru cci et alJ together with the visible data from 
DeMeo et a jj ( 2009I) . scaled to the visible albedo estimated by 



Lim et al. (2010). Overplotted are the basic model in red (#1 
in Table Q]) and model with additional 5% methane in blue 
(#2, ibid). Also shown are the visible and near-infrared col- 
ors of Vanth i n purp le. V, I and J band measurements are from 
iBrown et~aT ] (I2010I) . The colors of Vanth are normalized to 
Orcus' spectrum at 0.6 microns. 



First, to remove the dominant signature of the crystalline 
water ice from the spectrum we model the compos ition of Orcus 
using a code based on radiative transfer theory dHapkdfl993l) 
using optical constants of laboratory materials for inputs. We 
use optical constants of H2O ice i n both crystalline and amor- 
phous form (at 40 K and 38 K from lGrundv & Schmittll 19981 and 
Sch mitt et al.l 1 1 9981 resp ectively) and Titan and Triton Tholin 
dKhareet al.lll984l Il993h . The temperature of the H 2 optical 
constants are appropriate because the blackbody temperature at 
39 AU is about 43 K and Pluto's surface temperature ( with a sim- 
ilar s emi-major axis) is measured to be 40 ± 2 K dTrvka et al.l 
[1991 . Triton and Titan Tholins are used as representative ma- 
terial that aid in fitting the spectrum, because optical constants 
are available for these materials. However they could be replaced 
with different materials, such as other organics that have similar 
spectral properties. The models we present here (see Tab le Q] 
and F ig. [3) are based on the recent analysis by iDeMeo et alJ 
(2010) but differ slightly owing to the recent reevaluation of 
the albedo of Or cus from 20 + 0.03 dStansberrv et al.l 120081) 
to 0.27 ±0.06 bv iLim et all (I2010I) . We reduced the amount of 
amorphous water ice to about 10%, increasing the crystalline 
H20 by the same amount to better fit the 1 .65// b and and ad- 
justed the cosine asymmetry factor (Hapkel [l993l) to properly 
fit the data's higher albedo. The spectrum was then divided by 
the model 1, without methane (Table [TJ. We created a program 
in IDL designed to fit Gaussians to potential features in des- 
ignated wavelength regions. The program was set to search in 
the regio ns near 1.67, 1.72, and 2.2 pm, where methane absorbs 
strongly dOuirico &SchmitdlT997l) . We did not search for the 
band near 1.80 pm because of poor telluric correction in this 
wavelength range, nor the bands near 2.32 pm and 2.43 pm due 
to decreasing signal-to-noise ratio at wavelengths longer than 
^2.3 pm from low de tector sensitivity. A least-square minimiza- 
tion (IMarkwardtl2009l) was used to find the best-fit center, width, 
and depth of the bands. 

The results of the Gaussian fits are listed in Tabled and plots 
of the fits are shown in Fig. [4] We find a Gaussian fit near the 
1.67 pm feature, however, the center is at 1.654 pm indicating it 
is residual crystalline H2O not removed by the division of the 
data by the model. We do not find a band at 1.67 pm nor at 
1.72pm. The depth of the 2.209pm Gaussian fit is 9.5 ±2.3%. 



Table 1. Material present in the models of surface composition. 



Model 1 



Model 2 



Material 


Amount (%) 


Grain size 


Crystalline H 2 


60 


18 


Amorphous H 2 


10.5 


18 


Titan tholins 


2 


10 


Triton tholins 


6 


10 


Blue compound 


22.5 


10 


Crystalline H 2 


60 


18 


Amorphous H 2 


5.5 


18 


Titan tholins 


2 


10 


Triton tholins 


6 


10 


Blue compound 


22.5 


10 


Methane 


5 


100 



While we do not fit a Gaussian to the 2.32 pm feature we do note 
that the reflectance of the data is lower in this region with respect 
to the model that excludes CH4 (#1) suggesting there is a com- 
pound absorbing in this region, though not necessarily CH4. 

The question that remains is "Can the band at 2.2 pm be due 
to methane even though we do not detect other bands?" The rel- 
ative band depths of the 1.67 and 1.72 pm features compared 
to the 2.2pm feature are 95% and 88%, respectively, assuming 
a grain size of 200pm (the relative depths do not change sig- 
nificantly with variations in grain size of ± 100pm). We should 
clearly detect all three bands if the 2.2pm feature were due to 
methane only (as visible in Fig [3] by comparing the spectrum 
of Orcus with model #2, including 5% of methane to account 
for the 2.2pm feature), thus we place a limit of a maximum of 
about 2% of methane on the s urface. Our constrain t of 2% is 
lower than the constraint set in IDeMeo et alJ d2010h partly be- 
cause of the better quality of the data, but also because bands 
from 5% methane model at the higher albedo are weaker than 
for the lower albedo model. The majority (or all) of the 2.2pm 
absorption must therefore be due to either hydrated ammonia, 
ammonium, as suggested in previous work, or another yet un- 
known compound that absorbs in this region. Since H20:NH3 
and NH| do not have any other distinguishing features in the 
wavelength range of our data we cannot confirm or exclude their 
presence. 

We also search for bands at 2.27 4 and 2.314//m, the posi - 
tions of the strongest bands of ethane dOuirico & Schmitt» 1997). 
Ethane is a by-product of methane and its pote ntial presence 
on Or cus was suggested by models to data from Delsant i et al.l 
d2010l) . We do not find any bands at the locations of ethane ab- 
sorption within the limits of the quality of our data (which is 
the best to date). Either there is no ethane (less than about 5%, 
depending on grain size, within the limits of detection with our 
data) on Orcus' surface or, less likely, there is a concentration 
of ethane on the part of the surface observed by Delsan ti et al.l 
that we did not observe. We also check for the possible presence 
of other volatile compounds: CO2 (three narrow features near 
2pm), CO (1.578pm, we do not search at 2.352pm because 
of the poor quality of the data), N2 (2.15pm), and methanol 
(2.27 pm). We do not detect any of these features in our data. 



4.2. The colors of Vanth 

The spectrum of Vanth we extracted had a very low signal-to- 
noise ratio, therefore we drastically reduced its spectral resolu- 
tion and we only report in Table[3]its H- and K-band magnitude 
difference to Orcus. Fig.[3]compares the visible and near-infrared 
colors of Vanth to those of Orcus, scaled to the visible albedo of 
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Table 2. Parameters for weak bands. 



Species 


X e 


A m 


AA 


Depth 


10' - 




Qjm) 


Oum) 


(nm) 


(%) 




CH 4 


1.670 


1.654 ± 0.004 


8±5 


4.1 ± 1.9 




CH 4 


1.724 






< 2 




CH 4 


2.208 


2.209 ± 0.002 


6±2 


9.5 ± 2.3 




C 2 H 6 


2.274 






< 2 


5000 " 


C2H6 


2.314 






< 2 





Note s. B and centers (A e ) of the strongest CH 4 bands from 
lOuirico & Schmittl < ll997h . We list the band center (Am), width (AA), and 
depth, from the Gaussian fit to our data at these wavelengths if detected. 




52 1.64 1-66 1 
Wavelength (/j.m) 



1 .68 1 .70 1 .72 1 .74 1 .76 2.1 6 2.1 8 2.20 2.22 2.24 
Wavelength (/im) Wavelength (urn) 



Fig. 4. Gaussian fits to the data to search for weak bands. Two 
bands of approximately 4% and 9% are detected around 1.65 and 
2.2 //m respectively (see Table |2|. No band with depth weaker 
than about 2% is detectable around 1.72/mi. 

Table 3. Relative magnitude difference of Orcus and Vanth. 



B 



-5000 




-X [km] 

Fig. 5. Predicted position of Vanth relative to Orcus on 2010 
Fe b. 23.237, from the e xtrapolation of the initial 2004-2007 data 
bv lBrown et al . (2010). The two curves reprensent two possible 
solutions with symmetric pole solutions. The dots in the Eastern 
part of the frame correspond to the predicted positions at the 
time of the SINFONI observation for all retained orbital solu- 
tions given in Fig. [6] 



Filter 


Orcus 1 


Amag 


Flux ratio 


Refs. 


V 


19.36 ± 0.05 


2.61 ± 0.10 


9.0 ± 0.9 


1 


I 


18.63 ± 0.05 


2.31 ± 0.10 


8.6 ± 2.9 


1 


J 


20.64 ± 0.04 


2.66 ± 0.34 


11.5 ± 0.5 


1 


H 


21.65 ± 0.03 


2.34 ±0.12 


11.6 ± 1.3 


2 


Ks 




2.07 ± 0.24 


14.9 ± 3.2 


2 



References. (l) lBrown etaLl j201Cl) . (2) This work. 

Orcus from lLim et al.1 (|2010). Although the errors are large, the 
colors indicated that Vanth is both slightly redd e r in th e visible 
wavelength regime as reported by iBrown et al.l (1201 Oh and the 
near-infrared indicated by our data. While the colors of Vanth 
do not suggest a strong presence of water ice as seen on Orcus, 
the possibility cannot be ruled out within th e uncertainty of the 
data. As already addressed by IBrown et all it is not possible to 
draw any conclusions on the origin of this color difference to dis- 
tinguish between a capture or collision formation scenario. Even 
if Vanth is confirmed to have colors inconsistent with water ice 
it may not exclude a collision formation scenario. 

5. Orbital characterization 

The relative orbit o f the Orcus/Vanfh system has been deter- 
mined previously bv lBrown et al.l (12010 ) from observations con- 
ducted on the Hubble Space Telescope (HST), mainly with the 
High Resolution Channel of the Advanced Camera for Surveys 
(ACS), acquired over the period 2005-2007. Three years after, 
the uncertainty in the orbital period leads to a position uncer- 
tainty of Vanth along its orbit of about 15°, corresponding to 
2,000km (Pig. [5J. With one additional point for the relative po- 



sition of the components after 3 years, corresponding to some 
2,000 revolutions of the pair, it is possible to refine the orbit, 
especially the orbital period, which in turn allows refining the 
position prediction for future observations, in particular stellar 
occultations. 

We present in this section the relative astrometry measure- 
ments, the subsequent orbit computation, and the consequences 
on predictions of stellar occultations. 

5.1. Astrometry measurements 

To extract the relative astrometry of Orcus and Vanth, we first 
stacked the spectrocubes along the wavelength dimension to ob- 
tain two frames of higher SNR, one at each epoch. We then used 
the image model I presented in Sect. [3] the astrometry being 
given by the relative position of the Moffat functions centers. 

As visible in Fig.Q] Vanth is easily separable from Orcus on 
the image obtained in February, and we estimate the accuracy 
of its position to be about 25 mas, corresponding to a fourth of 
a pixel on the SINFONI detector. However, our observing date, 
which was imposed by the stellar appulse, was not optimal for 
taking full advantage of the parallactic effect and removing the 
ambiguity in the relative orbit inclination (see the two curves 
in Fig. [5). On the other hand, the moderate spatial resolution 
(FWHM of 380 mas) achieved in March did not allowed the lo- 
calization of Vanth, which flux was spread over a large area and 
diluted in the background noise (Fig. |T}. 

We also use the astrometric measurements given in 
IBrown et al.l (I2010I) . although, after an alyzing the public H ST 
data, it is clear that the dates given in IBrown et alJ d2010h are 
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wrong by a constant offset of half a day (confirmed by D. 
Ragozzine and M. Brown, personal communication). All the 
astrometric measurements are correctly reproduced in Table |4] 
along with the observation circumstances: heliocentric distance 
(A) and range to observer (r), right ascencion (a) and declination 
(5). 

5.2. Orbit improvement 

With the additional astrometric data we can improve the 
orbital parameters, reset the orbital phase of Vanth, and, tak- 
ing advantage of the time leverage, the orbital period is still 
better constrained. We computed an i mproved Keplerian orb it 
using a statistical inversion algorithm ( Hestroffer et al . 2005), 
which allows to probe a large portion of the orbital parameters 
space. We show in Fig. [6] all the solutions consistent with the 
observations, within the measurement uncertainties, and list the 
resulting improved orbital and physical parameters in Table [5] 

As mentioned before, the position of Vanth in its orbit 
was not favorable on 2010 Feb 23 (see Fig. [5]l to remove the 
ambiguity of the inclination of the orbit (+21 "with respect to the 
average plane-of-the-sky). However, the orbital pole solution 
that is almost perpendicular to the ecliptic (Pa, corresponding 
to i = +21°), coordinates ECJ2000 (321°, -2°), is statistically 
more probable given the observable data (see Fig. [5); although 
the other solution (Pb) inclined at about 52° from the ecliptic, 
coordinates ECJ2000 (340°, +38°), provides very similar 
residuals and cannot be ruled out yet. 

All the possible sol u tions w e find have low eccentricities 
(e<0.01). iBrown et af] (1201 Ol) had already highlighted the 
nearly circular nature of Vanth's orbit around Orcus, although 
we find here several solutions with larger eccentricities than 
the upper limit they reported (0.0036, 1 <x deviation to the 
data). Similarly, w e find an averag e orbital period l onger than 
IBrown et al.l ([2010) (consistent with O rtiz et aD201 ll however), 
their shorter period being still possible, although less likely. 
These differences can or iginate from the circular assumption 
on Vanth orbit used by IBrown et all In any case, since the 
eccentricity of the orbit is small, it suggests the tides have 
circularized the orbit, and subsequently that the satellite is also 
in synchronous rotation. The time for circularisation of such a 
system can be fast, of the order of r c ~ (a//?2) 5 « 1 ~ 10 6 years 
dFerraz-Mello et al.1 [2008). Another reasonable consequence 
from such tidal evolution is that the orbital plane is aligned with 
the equator of the primary. In this case, we would currently see 
Orcus under a small aspect angle (roughly 10° and 30° for Pa 
and Pb respectively). 

The circular orbit could also be due, independently of 

tidal effects , to a Kozai resonance ( Vl - e 2 cos i = constant, 
lKozaH[l962t) . in which case the inclination of Vanth's orbit with 
respect to Orcus' equator can be large. This means also that 
the direction of the spin for the primary is unconstrained and 
that its aspect angle can be much larger. We favor here the first 
solution (Pa) as it is slightly more probable statistically, and it 
is consistent with the low amplitude of the variat ions present in 
the o ptical lightcurve of Orcus (0.04 mag, see iThirouin et al.l 
I2010L and references therein). 

The total mass M is derived to high precision from Kepler's 
third law, considering this dynamical system is well isolated 
over the observations time scale, supported by the two-body 
fit. Then, t he rec ent size determination from ESA Herschel by 
iLim et all (l2010h allows us to evaluate the volume, hence den- 




0.005 0-01 0.015 0.02 0.005 0.01 0.015 0.02 

e e 

Fig. 6. Elliptical orbital elements: System mass (AT), semi- 
major axis (a), inclination (/), and period (P), plotted as function 
of the eccentricity (e). Each point represents a different orbital 
solution that fits the astrometry measurements in Table [4] See 
Table|5]for a summary of the orbital elements. 

sity, of the components. Because the components are not re- 
solved by the thermal infrared pixels, one can only derive a 
radiometric volume-equivalent diameter T> for the system, and 
subsequently model-dependent values for the density. This ef- 
fective diameter can be simply related to the volume-equivalent 
diameter of each component by considering their projected ap- 
parent surface D 2 -D 2 +D 2 . Assuming that the two components 
are of similar density, one eventually gets a rough estimation of 
the bulk density: 

_ 6_M (1 +7? 2 ) 3/2 
9 ~ (1 +5R 3 ) ( } 

where the size ratio of the bodies % = DJD a = 1 o-°- 2A '«, 

essentially depend on their apparent magnitude difference 
(Am » 2.4 ± 0.4, see Ta ble [5), modulated by their unknown 
albedo ratio. Following Bro wn et al ] (I2010h . we consider two 
possible albedos for Vanth: 0.27 (i.e., equals to Orcus', corre- 
sponding to co-formation mechanisms, Section l5~3l l. and 0.12, 
typical for TNOs exempt of water ice (for capture mecha- 
nisms, Section 15.3b . We find mass ratios of 30 and 8, for 
the equal-albedo and diffe rent-albedo a ssumptions re spectively, 
close to the estimates by IBrown et alJ d2010h and lOrtiz et all 
d201lh . Based on the re fined volume-e quivalent diameter of the 
Orcus/Vanth system bvlLim et al.l (120101) . which is 10% smaller 
than the estimate by IBrown et al.l (1201 Ol) . we find a density of 
about 2.3 ± 0.8 g/cc. In any case, the uncertainty on the bulk 
density is dominated by the uncertainty on the size determina- 
tion, which cannot be more accurate than about 10% dLim et al.l 
2010). The relative precision on the bulk density cannot there- 
fore be smaller than about 30%. Better knowledge of the spin 
properties of Orcus, and of its shape are now required to improve 
its size, hence volume, estimate. 

5.3. Formation Mechanism 

Trans-Neptunian binaries can be used to impose con- 
straints on various models of Solar System formation 
evolution, e.g., to constrain N eptune's migration history 
dMurrav-Clav & Schlichtingl 1201 ll) . Several formation mecha- 
nisms have been proposed and each leads to different predictions 
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Table 4. Astrometric data for the relative position of Orcus/Vanth used to reconstruct their mutual orbit. 



Date 1 ' 


AX* 


A* 4 


A 






s 


Instr. 


(JD) 


(") 


(") 


(AU) 


(AU) 


(°) 


(°) 




2453687.66400 


0.206 


-0.147 


47.811 


47.702 


144.039 


-4.356 


ACS 


Z404U4U. joyUU 


U.ZZD 


n 1 1 1 
-U. Ill 


A Q fY7< 
4o. U/ J 


Al 1^1 


1/1/1 Q 1 ^1 


A TOO 
-4. 1 Ay 




2454044.36600 


-0.258 


-0.005 


48.013 


47.752 


144.839 


-4.766 


ACS 


2454051.57900 


-0.006 


-0.243 


47.898 


47.753 


144.871 


-4.828 


ACS 


2454056.08900 


-0.036 


0.240 


47.824 


47.754 


144.882 


-4.865 


ACS 


2454066.14600 


0.053 


0.240 


47.659 


47.755 


144.882 


-4.938 


ACS 


2454080.33600 


-0.030 


-0.244 


47.433 


47.757 


144.824 


-5.019 


ACS 


2454416.29300 


-0.263 


-0.024 


47.968 


47.801 


145.703 


-5.304 


NICMOS3 


2454439.78000 


0.245 


0.078 


47.585 


47.804 


145.696 


-5.473 


WFPC2 


2455250.73686 


0.260 


0.052 


46.954 


47.895 


146.287 


-6.396 


SINFONI 



The first block is reproduced from lBrown et al.l |2010) with the corrected values for the dates (see text). The second bloc correspond to our 
appulse observation on 2010 February 23. 
Positions are positive through East and North. 



Table 5. Orbital elements and physical parameters for the 
Orcus/Vanth system. Uncertainties are the formal 1 <x deviation 
(from Fig. [6]l and do not account for possible systematics (that 
arise from model incompleteness, like, e.g., neglecting Orcus J2 
which is unknown). 



Qty. 


Value 


Uncert. 


Description 




a 


9 030 


89 


km 


Semi-major axis 


e 


0.007 


0.003 




Eccentricity 


i 


(±)21 


2 


deg. 


Inclination 1 


Q 


144 


63 


deg. 


Ascending node 1 


n 


37.732 


0.002 


deg/day 


Mean motion 


P 


9.540 6 


0.0004 


day 


Orbital period 


Pa 


(321,-2) 


3 


deg. 


Pole A coordinates 2 


Pb 


(340, +38) 


3 


deg. 


Pole B coordinates 2 


M 


641 


19 


10 18 kg 


Total mass 


O 


850 


90 


km 


Diameter 3 


D 


807 or 761 


100 


km 


Orcus diameter 4 


D, 


267 or 378 


100 


km 


Vanth diameter 4 


P 


2250 or 2470 


885 


kg irr 3 


Bulk density 



1 In the tangent plane. 

2 Pole coordinates are given in the ecliptic J200 reference frame. 

3 Volume-equivalent from radiometry (see lLim e t al. 2010). 

4 Assuming alternatively p v = p a , or p = 0.27 and p v = 0.12. 



for the physical and orbital properties of TNBs. These forma- 
tion mech anisms can, broadly, be broken down into fo ur classes: 
collis ion dWeidensc hilling 2 0021). rotationa l fission (|Ortiz et al.l 
201 ll), capture (e.g.. IGoldreich et ail 120021: iFunato et alj 120041 
Asta khov etail 120051: iLee et al] 120071: iGamboa Suarez et al] 
2010), and gravitational coll apse dNesvornv et al.ll2010h . 

In the collision model of lWeidenschi lling (2002) two objects 
collide inside the Hill sphere of a third object. These objects then 
fuse into a single object thereby p roducing a binary. The capture 
models of IGoldreich et al] d2002l) rely on two objects interpen- 
etrating their mutual Hill sphere and then being stabilized ei- 
ther through dynamical friction (the L 2 mechanism) or through 
a scattering event with a t hird, s imilarly sized, object (the L 3 
mechanism). IFunato et al.l (|2004) proposed a hybrid collision- 
capture mechanism. Initially two objects collide to produce a 
binary whose components (as is usual for a collision) have quite 
different masses. Subsequently, exchange "reactions" with larger 



third bodies displace the secondary and so ramp up the mass ra- 
tio. This eventually leads to binaries having similarly sized part- 
ners. However, this mechanism appears to lead to orbital prop- 
erties (in particul ar, ellipticities) dissimilar t o those actually ob- 
served (e.g.. lNollll2003l:lAstakhov et alj|2005l) . 

A further capture scenario, chaos-a ssisted capture (CAC, see 
Astakhov et al. I l2005t lLeeeTSll2007b . was originally proposed 
to ex plain the capture of i rregular moons at the giant planets 
(e.g JAstakhov et al.ll2003l) . The scenario for the formation of 
TNBs is as follows: two objects initially become caught up in 
very long living, yet ultimately unstable, chaotic orbits within 
their mutual Hill sphere. During this phase the binary may be 
permanently captured and subsequently hardened through mul- 
tiple scattering encounters with relatively small "intruder" bod- 
ies. The CAC model predicts similarly-sized binary components 
whose mutual orbits are in good agreement with observations. 
In fact, it was proposed that the propensity to form similar sized 
binary partners is a direct fingerprint of chaos assisted capture. 
A further specific prediction of the CAC mechanism concerns 
the distribution of retrog rade to prograde mutual orbits, i.e., mu - 
tual orbit inclinations (Asta khov et al] 120051: ILee et al] 120071) . 
In particular, retrograde mutual orbits are predicted to be rel- 
atively common. The C AC mechanism has been criticized by 
ISchlichting & Saril d2008l) who argue, e.g., that formation via 
transient, chaotic binaries is not as important as the L 2 and L 3 
mechanisms. These objectio ns are addressed in some detail by 
Gam boa Suarez et al. (20lj. 

Finally. iNesvornv et al] d2010h recently proposed a model in 
which TNBs formed during gravitational collapse. Angular mo- 
mentum considerations in the planetesimal disk explain the for- 
mation of binaries rather than condensation into a single object. 
This gravitational instability model predicts identical composi- 
tions and colors for TNB partners and also inclinations generally 
i < 50°, i.e., retrograde mutual orbits are predicted to be rare. 
Unfortunately true inclinations of mutual T NB orbits is currentl y 
available for a handful of systems only (see Gru ndv et al.ll201 ll) . 
For example, Table [5] and Figs. [5] and [6] show that only knowl- 
edge of the true inclination of the mutual orbit (i.e., by making 
further observations to remove the symmetry ambiguity) will re- 
veal if the mutual orbit of the Orcus/Vanth system is prograde or 
retrograde. 

We remark that the CAC model also allows for the direct for- 
mation of almost circular orbits similar to that of Orcus/Vanth. 
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This mechanism was proposed as a possible origin of the appar- 
ently almost circular orbit of the T NB 2001 QW322 dPetit et al.1 
2008; iGamboa Suarez et al.l 120101) . It is conceivable that the 
Orcus/Vanth system was produced in a similar manner. In sum- 
mary, determining the colors and compositi ons, as well as the 
mutu al orbital true inclinations of TNBs JSch lichting & Saril 
2008) remain important quests. Knowledge of these properties 
will provide important information to constrain, and to distin- 
guish between, the variously proposed TNB formation models. 

5.4. Stellar occultation 

Observations of stellar occultation phenomena are of high inter- 
est to derive, in a direct and most precise way, the size of the 
occulting body. When several chor ds are gathered one can de- 
rive the projected shape contour (see lMillis & Dun ham 1982, for 
a extensive summary). Moreover such observations can reveal 
presence of a teneous atmosphe re, and then derive its pressure 
profile (e.^. lSicardvet al.ll2003l) . 

In the case of a binary system it is also important to be able 
to predict the pat h of each compone nt more than the one of the 
centre of gravity (Assafi n et al.l20l6T) . Given the large separation 
of Vanth with respect to Orcus (~9 000 km), the prediction of an 
event by the secondary can be either deplaced in time by sev- 
eral minutes when the binary system is oriented in the direction 
of apparent motion, or fully decoupled from the one of the pri- 
mary otherwise. Indeed, the occultation path by the primary may 
not be visible on Earth while the one by the secondary could be 
visible. 

Owing to the current uncertainty on TNOs ephemerides, pre- 
cise prediction of occultation requires "last minute" astrometry 
measurements with respect to background star to obtain reliable 
predictions of the occultation track on Earth. Besides, the rela- 
tive position of the secondary is needed to a precision of approxi- 
mately 1 00 km to enable putting in place useful observation cam- 
paigns. Given the improvement on the period estimate, and the 
orbital phase reset provided by our astrometric point obtained in 
2010, we can derive better prediction for the location of the sec- 
ondary and subsequently better prediction for stellar occultations 
by Vanth. Given the possible mass-ratio range the path on the 
Earth of the primary can also be shifted from the position of the 
centre-of-mass. Assuming the components have same albedo the 
mass ratio is of the order of «0.02-oce, corresponding to a shift 
of 180-iso km. We list in Table [6] the relative position of Vanth 
around Orcus for upcoming stellar occultations. The relative po- 
sitions in Table [6] corresponds to the nominal values; the asso- 
ciated error distribution is not exactly Gaussian, but an ellipse 
curved and stretched along the trajectory (similar to Fig. [5j. The 
one sigma uncertainty is large, approximately ±450 km along 
the secondary's nominal trajectory (about +120 km across), for 
both occultations. Depending on the occultation occurence (mo- 
tion relative to the star and orientation of this uncertainty ellipse) 
it can correspond to a shift of the occultation path on Earth (or 
miss) or a shift in time. Updates are available onlinfl 

6. Conclusion 

From spectro-imaging measurements of the Orcus-Vanth system 
we were able to obtain a spectrum of Orcus, H and K magnitude 
for Vanth, and determine their mutual position. We searched for 
possible weak bands in the spectrum of Orcus, limiting the pos- 
sible amount of CH4 to no more than ~2%, and ethane to ~5%. 

2 http://www.lesia.obspm.fr/perso/bruno-sicardy/ 



Table 6. Relative positions (6x, 6y) of Vanth with respect to 
Orcus for future occultation events (6x positive through East). 
Uncertainty on position draws an ellipse of ±450x120 on 
Earth, around the reported positions (see text for details). 



Date 


<5x (km) 


Sy (km) 


2012 Mar 21 


-6977 


-5259 


2012 Nov 16 


-8937 


347 



The presence of other compound(s) is required to explain the 
depth of the absorption band detected at 2.2 fim. Presence of hy- 
drated ammonia, or ammonium, could explain this depth, but the 
lack of distinguishable features for these species in the observed 
wavelength range forbid any strong conclusion to be drawn. 

Vanth appears slightly redder than Orcus in the visible to 
near-infrared wavelength range, although uncertainties in the 
measurements are large. While this suggests that Vanth's com- 
position, or age, is different from Orcus, it does not constrain 
formation scenarios (collision, co-formation, or capture). By de- 
tecting Vanth a t a new posi t ion th ree years after previous orbit 
calculations bv lBrown et al. I J20101) we are able to reset the phase 
of the orbit, allowing prediction of Vanth's track on Earth during 
future stellar occultation events. 
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